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ABSTRACT 

This paper investigates the roles of pressure-strain and turbulent diffusion models in the 
numerical calculation of turbulent plane channel flows with second-moment closure models. 
Three turbulent diffusion and five pressure-strain models are utilized in the computations. 
The main characteristics of the mean flow and the turbulent fields are compared against 
experimental data. All the features of the mean flow are correctly predicted by all but one 
of the Reynolds stress closure models. The Reynolds stress anisotropies in the log layer are 
predicted to varying degrees of accuracy (good to fair) by the models. None of the models 
could predict correctly the extent of relaxation towards isotropy in the wake region near 
the center of the channel. Results from the direct numerical simulation are used to further 
clarify this behaviour of the models. 


*This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Second-moment turbulence closure models first appeared about four decades ago with the 
proposal of a simple linear model for the pressure-strain correlation by Rotta (1951). Models 
with increasing complexity and sophistication have followed. The milestones are the works of 
Daly and Harlow (1970) (denoted DH), Hanjalic and Launder (1972) (denoted HL), Launder, 
Reece and Rodi (1975) (denoted LRR), Shih and Lumley (1985) (denoted SL), Fu, Launder 
and Tselepidakis (1987) (denoted FLT), and Speziale, Sarkar and Gatski (1991) (denoted SSG). 
Speziale (1991) has presented an insightful review of this hierarchy of closure models. The 
earlier models (DH, HL and LRR) propose approximations for the pressure-strain correlations 
which are linear in the Reynolds stresses, whereas the latter models (SL, FLT and SSG) use non- 
linear terms for these correlations. Although second-moment closure models have a sounder 
theoretical basis than eddy-viscosity based models such as the k-e model, their superiority in 
predictive ability has not been demonstrated in a consistent and systematic manner. Most of the 
models were derived with the use of homogeneous flow assumptions and the original applications 
have emphasized the prediction of homogeneous or nearly-homogeneous flows. There is a need 
for systematic studies in which the models are applied to flows with increasing complexity 
which are of practical importance, and the computed results compared to experimental data and 
results obtained with two-equation models, with the goal of establishing their capabilities and 
inadequacies. It is more usual to see a quantum leap in the application of the models to calculate 
highly complex two- and three-dimensional flows (Amano and Goel, 1987; Sykes et al., 1986; 
Demuren, 1992). From the point of view of developing turbulence models, such an exercise is 
usually inconclusive since the performance of the models in the (inhomogeneous) elementary or 
component flows is not well known. 

The present study is the first stage of an attempt to bridge this gap. The pressure-strain 
models tested here are : the quasi-isotropic model of LRR with and without wall-reflection 
terms, the models of SSG, SL and FLT, all without any special wall-proximity treatment. Three 
formulations for the diffusion terms are also examined. These are proposals by DH, HL and MH 
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(Mellor and Herring, 1973). The first has been criticized for violating the symmetry of indices 
in uiujuk (the triple velocity correlation) but is still widely used because of its simplicity. The 
last two preserve this symmetry in the indices of uiujuk. 

The test problem is the fully-developed plane channel flow at high Reynolds number. 
Surprisingly, it was difficult to find complete sets of experimental data which fulfill the following 
requirements ; high Reynolds number, high aspect ratio, long development length and high 
accuracy and consistency. Comte-Bellot (1965) presented data at very high Reynolds numbers 
but these showed an inexplicable Reynolds number dependence. Channel flow data by Clark 
(1968) also showed an excessive dependence on Reynolds number, and values of the normalized 
turbulent kinetic energy k/U£ appeared to be too high in the near wall region. Laufer (1951) 
presented data for a 12:1 aspect ratio channel at Reynolds numbers (based on bulk mean velocity 
and half- width) in the range 10,300 to 52,000 taken at 86 half- widths from the inlet. Hussain 
and Reynolds (1975) suggested that this length would be insufficient for full development of the 
turbulence field. They presented measurements at 450 half-widths to support this, but their data 
set was incomplete since only the longitudinal normal stress component was given. A review 
of duct flow measurements by Klein (1981) supports their assertion. Nevertheless, Laufer’ s 
data represent then the best compromise, and hence was used to construct the anisotropy stress 
tensor components to which the present model computations are compared. It was supplemented 
with direct simulation results of Kim, Moin and Moser (1987) (denoted as KMM), and Kim 
(unpublished), in order to explore more detailed features of the turbulence fields in comparison 
to the turbulence models. 

2. Mathematical Formulation 
2.1 Mean Flow Equations 

The Reynolds-averaged mean-flow equations for steady, incompressible turbulent flow can 
be written in Cartesian tensor notation as: 
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Continuity 


( 1 ) 


dUj 

0Xj 


= 0 


Momentum 


a lv ,,, iap, A 

W ' «x, 


n/dUi du,\ 

h iU, + + 


( 2 ) 


where x;=(xi, X 2 , X 3 ) represent the Cartesian coordinates, Ui=(Ui,U 2 ,U 3 ) the Cartesian mean 
velocity components, P is the pressure, p the molecular viscosity and p the density. Einstein s 
summation rule for repeated indices is utilized. u;u/ (with i=l,2,3, and 1 = 1 , 2,3) represents the six 
components of the Reynolds stress tensor, Rjj which must be determined by the turbulence model. 


If, in the fully-developed plane channel flow, the walls are in the (1 — 3) planes as shown in 
Fig. 1, then, ^ = 0 in equations (1) and (2). Also, the Reynolds stress tensor R a will 

have only 4 non- zero components; the three diagonal elements u^,u^ and uj which represent 
the normal stresses, and one off-diagonal element ujul which represents the shear stress. 


2.2 Reynolds Stress Equations 

The transport equations for the Reynolds stress components can be written for high Reynolds 
number turbulent flow in Cartesian tensor notation as: 

(Ufujuj) = Dij + P ij + 7Tij — fq (3) 

where Dy is the diffusion, Py is the production, 7ry is the pressure-strain correlation, and ey is 
the dissipation rate. 

The production term is Py = -u[u - u]u 7 and the dissipation is assumed to be 
locally isotropic so that ey = 2/3 6 tJ e, where e is the dissipation rate of the turbulent kinetic 
energy k, to be determined from the solution of a transport equation, and 6 y is the Kronecker 
delta. However, FLT attempts to account for anisotropy in the dissipation, as will be discussed 
in a later section. 

The models for Dy and 7 ry are the subjects of the present study. 
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2.3 Diffusion Models 


The three diffusion models examined are based on the proposals of DH, HL, and MH. They 
are written in tensor form as, 


Dii = (-%) 




(4) 


where ( ),k represents the first derivative with respect to Xk, and Ty k is given by: 


DH 


: -Tijk = c s i ^ J^ui t u/(uiuj) )1 


(5) 


HL : 


-T 



iU^UjUk)^ 4- UjU,(u k U,) , + U k U/(UiUj) , 


( 6 ) 


MH : — Tiik = 


k 2 


1C 

ijk = c s3— ( u i u j) k + ( u i u k)j + ( u j u k) } 


(7) 


with c s i=0.22, c S 2=0.11 and c s 3=2/3c s2 . These are essentially gradient diffusion models in which 
DH and HL have non-isotropic diffusion coefficients, but that for MH is isotropic. HL and MH 
diffusion models preserve symmetry in the indices but the DH model does not. In a general 
three-dimensional flow, the HL model requires the evaluation of 27 derivatives, whereas the DH 
and MH models each contain only 9 derivatives. For developed thin shear flow with ( ) a = 
( ), 3 =0, the diffusion terms are assembled in Table 1. These models result in expressions for 

the diffusion of k which are different from those usually given in two-equation models. 


2.4 Pressure-Strain Models 


Five models for the pressure-strain correlation are examined in this study. These are 
the quasi-isotropic version (model 1) of LRR, with and without wall-proximity treatment, the 
dynamical-systems based model of SSG, and the non-linear models of SL and FLT, both of which 
are derived from realizability constraints. The LRR model without wall treatment is denoted by 
LRRNW in this paper. In their model 1 , LRR proposed to account for wall-proximity effects 
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by making coefficients in the equations functions of the average distance from walls. In the 
simpler version (model 2) usually called the Gibson-Launder model, they are treated as wall- 
reflection terms. It should be noted that these wall-reflection terms are not near-wall corrections 
in the conventional sense, since they are applicable to the fully turbulent region beyond the 
viscous sublayer and the buffer zone, and they still have significant contributions at the center 
of the channel. However, there is uncertainty as to how rapidly the functions should decay with 
distance from walls, or how to estimate average distance from walls in complex geometries. 
Hence, it is now generally accepted that the need for wall-proximity treatment is an undesirable 
feature in a pressure-strain model. The comparison between LRR and LRRNW serves merely 
to show the dominant role wall reflection terms may play in determining the anisotropy levels. 

The primary question that we address in this study is how do the sophisticated Reynolds 
stress models perform in a ‘building block’ inhomogeneous shear flow such as channel flow. 
We are interested in evaluating the performance of these high Reynolds number models in a 
region (y+ > 200) away from the wall - an issue separate from near-wall turbulence modeling. 
Therefore, the same wall-function treatment is used to bridge the near-wall region in all the 
models considered here. 

The pressure-strain models can be written in terms of the anisotropy tensor 
by (= ufuj/2k — jtfjj) , the rate of strain tensor S,j , the rotation tensor 

Wij (= 1^- - , and the rate of production of turbulent kinetic energy P k in the general 

form : 


TTij — aoebij -|- ot\c(^b xk bj k l/3II<5jj) T cx2^^ij + + 

+ k{ a 4 {bikSjk + b jk S tk — 2/Z6,jb kt S kl ) + ot${b xk W jk -f b jk W lk ) 


linear terms 


+ a 6 {bikbikSji + b ]k bi k S,i - 2b kj bi t S k i) + a7(bi k b lk Wji + bj k bi k W t i) 

V -v- y 

quadratic terms 

\pnl^ln {pik ^ j k bjfcW » A: ) "I" ^bniibnj^mk^^nk ] } 


cubic terms 
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The model coefficients <*o . . . a& may be, in general, functions of the invariants of the anisotropy 
tensor. The corresponding relations for the five pressure-strain models are presented in Table 2. 

In the table, f is a wall-proximity function which takes a value of unity in the fully turbulent 
region near a wall and zero in a flow free from walls. LRR proposed a linear decay for f, 
but Demuren and Rodi (1984) found the wall effect too strong near the center of the channel, 
hence they prefer the quadratic form used in this paper. II(= b^by) and III(= bikbic m b m i) 
are, respectively the second and third invariants of the Reynolds stress anisotropy tensor and 
Rt is the turbulent Reynolds number. The constants c^ and k take the standard values 0.09 
and 0.42, respectively. 

The first line in equation (8) contains a mixture of terms representing both “slow” and “rapid” 
contributions to the pressure-strain correlation. The first term is the usual Rotta term for the return 
to isotropy. All the models have this term, which has a constant coefficient in the LRRNW and 
SSG models, but is a function of the wall-proximity variable in the LRR model, and a function of 
the invariants in the SL and FLT models. In the latter, the particular forms are derived to ensure 
that the turbulence remains realizable in the two-component limit as the wall is approached. 
Such a condition is never approached in the present study since the integration is not performed 
all the way down to the wall. In the FLT model, unlike all the other models, this term does not 
contain the full return to isotropy term, since the value of the coefficient o-o may become less 
than 2 when the second invariant II gets very small. The remaining contribution, which ensures 
the return to isotropy is contained in their approximation for e tJ 2/3 + 

The treatment corresponds more nearly to the usual practice if the second part of this expression 
is combined with the return to isotropy term. Only the SSG and FLT models have a non-linear 
contribution to the return to isotropy. The last two terms on the first line are contributions to the 
“rapid” part, the first of these is a linear term and the second is quadratic in by, since P* itself 
is linear in by. The major contributions to the “rapid” part of the pressure- strain correlation are 
those collection of terms in lines 2 to 4 of equation (8). All models contain the linear terms. 


2 _ 2F ly/2 


bijA . 
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Only the SL and FLT terms contain the quadratic terms and only the FLT model has terms 
which are cubic in the anisotropy stress tensor. 

2.5 k-e Model 

Calculations were also made with the standard high-Reynolds number form of the k-e 
turbulence model. The equations for k and e can be expressed in tensor notation as: 

^-(U,k) = D k 4- P k - e (9) 


d 2 

^-(U,e) = D f + Q, ^P k - c (2 ~~ (10) 

In the standard form of the model the terms D k and D e are approximated by gradient diffusion 
relations as: 


D k = 



( 11 ) 


D e = 



( 12 ) 


These equations are routinely solved, even when the interest is only in the solution with the 
second-moment closure models, e is of course required for closure. The trace of the Reynolds 
stresses should be equal to 2k, so that the solution of the k-equation is redundant. It serves, in 
this study, solely as an additional check for the convergence and the consistency of the solution. 
Thus, it was required that half the trace of the computed Reynolds stress distributions should 
agree with the computed distributions of k to within 0.2%, at every point in the flow. For 
consistency in this case, D k must take the forms given in Table 1. The emphasis in the present 
study is on the models for the pressure-strain correlation and the turbulent diffusion, so the use 
of anisotropy diffusion coefficients in the e-equation, as proposed by HL may becloud the issue. 
The empirical constants in equations (9-12) are : <r k = 1.0 ; a c = 1.3 ; c c i = 1.44 ; c c 2 =192. 
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Note that some authors (for example SL) use versions of the k-e equation with different 
values for the model coefficients. The question arises as to whether such differences change the 
equilibrium Reynolds stress anisotropies in the log-layer. Consider the Reynolds stress equation 
in the fully-developed, incompressible channel flow 

Dij + Pij + TTij — g c ^ij = 0 (^3) 

In the log-layer, where diffusion can be neglected and (13) becomes 

2 

P ij + - 3 e ^i = 0 

or 

= o (14) 

t e 3 

All the models for Try are of the form = /(by, in a simple shear flow. Since P k = e in 
the log-layer, we have ^ , and thus = /(by). Similarly, ^ = g{b y, §*) = Sf(by). 

Thus, equation (14) becomes an algebraic equation for by independent of e. Consequently, the 
anisotropy tensor by in the log- layer is independent of the form and coefficients of the e ^equation. 
Note that by is determined completely by the form of the pressure-strain model and should be 
constant if there is no explicit y dependence in the model. 

2.6 Solution Procedure 

The problem, as formulated, is strictly one-dimensional, but a two-dimensional TEACH-type 
code which solves the full, time- averaged, Navier-Stokes equations is utilized, in order to ease 
extension to other flow cases in future studies. The redundant terms are simply set to zero in the 
present study. The initial conditions for the mean flow and turbulent stresses are taken from the 
simulation data of Kim, Moin and Moser (1987), but these are scaled up to yield an effective 
Reynolds number (based on bulk mean velocity and half width) of 5.2x10^, to coincide with 
the highest Reynolds number of Laufer’s experiments. The computed results (normalized with 
U r ) are similar to those obtained (Demuren and Sarkar, 1991) at a higher Reynolds number (~ 
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4 x 10 5 ). Computations are performed for half of the channel from the lower wall to the mid- 
plane. The calculations use 32 grid points in the transverse direction and sufficient lengths (over 
400 half-widths) in the longitudinal direction to ensure full flow development. Computations 
with twice as many points in each direction do not produce significantly different results. The 
objective of this work is to compare the behavior of various pressure-strain and diffusion models, 
so the computations are for the high Reynolds number flow region only in which the viscous 
sublayer is not resolved but is bridged using the standard wall-function method. Along the line 
of nodes nearest to the walls (y + ~30) local equilibrium is assumed: The streamwise velocity 
component is specified based on the logarithmic velocity of the wall (U + = £ In y + + 5.0) ; k 
= U r 2 / Cfl W ; e = U T 3 /(«y) ; = 1.07 k ; TIJ = 0.41 k ; = 0.52 k ; = - 0.30 k. The 

specification of the second set of boundary conditions is that the first derivative of all dependent 
variables is set to zero normal to the symmetry plane. 

3. Results and Discussion 

The results of the computations to test the turbulent diffusion models are compared to results 
obtained with the k-e model and experimental data in Figs. 2—4. For each of the computations 
with the second-moment closure the pressure-strain correlations are modelled with the SSG form, 
but the three diffusion models given by equations (5-7) are used. Profiles of the streamwise 
velocity are presented in semi-logarithmic form in Fig. 2. All the profiles agree very well with 
the log-law (U + = £ In y + + 5.0) in the inner layer. The results with the k-e model display a 
pronounced wake region, but the second-moment closures give only a small wake region. Since 
the channel flow has a favorable pressure gradient, the wake region is of course much smaller 
than that for a flat plate boundary layer. Measurements by Clark (1968), Laufer (1951) and 
Hussain and Reynolds (1975) also show a very small wake region. Hence, the velocity profiles 
are well predicted by all the diffusion models. Another test of the diffusion models is in the 
prediction of the Reynolds stress anisotropy tensor in the relaxation region 0.7> y/6 < 1.0. Figure 
3 compares these results with the data of Laufer for the 4 non- zero components. As seen in the 
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1 1 and 22 components, the MH model produces the strongest rate of relaxation towards isotropy 
as the center of the channel is approached, and it gives the best agreement with experimental 
data. The DH model shows little change in the anisotropy level between the log-layer and the 
center of the channel. The results for the HL model lie somewhat between the other two. For 
the 33 component, the DH model shows slightly better agreement with the data. These results 
can be explained by considering order of magnitude estimates of the diffusion terms given in 
Table 1. DH and HL models have coefficients c s i and c S 2 which are respectively 3 and 1.5 times 
larger than that for the MH model. But they are modified by anisotropy coefficients, the most 
significant of which is u\fk, with a value 0.4 in the log-layer. Order of magnitude estimates 
then show that the diffusion terms for uj are roughly equal, but for u\, the MH model yields a 
diffusion term which is about 2.5 times greater than that for the DH model and 1.7 times greater 
than that for the HL model. The budget of the Reynolds stresses from the present study (not 
shown) and from the DNS data of Kim (unpublished), shows that the diffusion of the normal 
stresses are positive in the central part of the channel (y /6 > 0.5) while the diffusion of the shear 
stress is negative. Thus, an increase in the diffusion of u\ would increase an d move b 22 
closer to zero, thereby reducing the anisotropy, if uj and u\ are not correspondingly increased. 
This is the case here, since the diffusion models yield magnitudes which are more nearly equal 
for u\ and u\. Therefore, the MH model which has the largest diffusion term for produces 
the fastest rate of relaxation towards isotropy and the best agreement with experimental data. 
Similar order of magnitude estimates can be used to explain the results for b 33 and b^. Since 
the diffusion of uiui is negative, a higher magnitude will produce a lower value of bi 2 . The 
comparisons of the predicted turbulent kinetic energy are presented in Fig. 4. The predictions 
are in reasonably good agreement with the data. The slight differences between the predicted 
results can also be explained by the aforementioned order of magnitude estimates. 

The predictions with the five pressure-strain models are compared in Figs. 5-7. The MH 
model for the turbulent diffusion is used in each case. Figure 5 shows that all models except 
the SL model give reasonable prediction of the mean streamwise velocity profile in agreement 
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with the universal logarithmic law of the wall. The explanation for this can be found in Fig. 
6 and Table 3. The main requirement for the correct prediction of the mean velocity is that 
the model should yield an accurate distribution of the shear stress. The latter is related to the 
bi 2 component of the Reynolds stress anisotropy tensor. Table 3 shows that this component is 
reproduced fairly accurately by all but the SL model. In the inner core of the flow (y /6 < 0.8) 
the SL model underpredicts bi 2 by about 30%. Now, in the equilibrium layer for thin shear 


flows the ratio of production to dissipation is given by : 



(15) 


and since P k/e is approximately unity the normalized shear rate (Sk/c) is inversely proportional to 
bi 2 . Underprediction of bi 2 will produce excessive shear rate and hence a poor velocity profile. 
Experimental and DNS data suggest that Sk/e should be equal to about 3.3 in the inner layer. 
Most of the models predict values in the range 3.0 — 3.5, but the SL model predicts values 
of around 4.3. The latter is not surprising since Speziale and Mac Giolla Mhuiris (1989) had 
reported that the SL model predicts equilibrium values, in homogeneous shear flow, of bi 2 and 
Sk/e of -0.12 and 6.93, which are respectively, lower and higher in magnitudes than experimental 
values by about 20%. The channel flow is of course not a homogeneous shear flow but there 
are some similarities. For example bij and Sk/e have constant values in both the log-layer of the 
channel and the equilibrium homogeneous shear flow. 


The comparisons of the normal components of the Reynolds stress anisotropy tensor in Fig. 
6 are even more instructive. A comparison of results obtained with the LRR model to those with 
the LRRNW model shows the effects of the wall-reflection terms. They produce a significant 


increase in the anisotropy of 11 and 22 components, with little effect on the 33 component. 
Although the increase in anisotropy is strongest near the wall, it remains pronounced even at 
the center of the channel. This is contrary to expectation. The dilemma is how to devise a 
function for the wall reflection effects which decays at the right rate away from the wall that 


would also be general enough for application to more complex flows. Better still, the model 


should not require wall reflection terms. 
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If we consider the level of anisotropy in the inner layer, the SSG model (dash-dot curves 
in Fig. 6 ) gives the closest predictions of bn, bn and 1 ) 33 . The LRR model gives reasonable 
predictions for bn and b 22 . but not for b 33 . Both these models are relatively simple, without 
the quadratic and cubic terms in equation ( 8 ). Their good performance is probably due to their 
superior calibration for homogeneous flows. The SSG model was calibrated in a dynamical 
systems approach using data from the nearly-homogeneous shear flow experiments of Tavoularis 
and Corrsin (1981). The LRR model was calibrated with the earlier (and probably less reliable) 
experimental data for nearly-homogeneous shear flow of Champagne, Harris and Corrsin (1970) 
and the wall-reflection part of the model used a consensus of near-wall data. It is surprising 
that the FLT model which contains both the quadratic and cubic terms in equation ( 8 ) does 
not give predictions which are superior to the much simpler models. Furthermore, b 33 shows 
an increase in anisotropy towards the center of the channel. Perhaps the calibration of the 
model is to blame. Numerical experiments indicate that the level of anisotropy in the 1 1 and 
22 components can be increased by increasing the magnitude of the coefficient of the cubic 
terms, as in equation ( 8 ). However, this produces little effect on the 33 component. Again 
the predictions by the SL model of bn and b 22 in the channel flow are much lower than in the 
experiments. Correspondingly, b 33 is much too high. Such underprediction of bn and b 22 has 
also been observed in homogeneous shear flow (Speziale and Mac Giolla Mhuiris, 1989). The 
values of the Reynolds stress anisotropy tensor in the log-layer are summarized for all models 
and the experimental data in Table 3. 

A notable feature of the predictions is that the models do not fully reproduce the rapid 
relaxation towards isotropy (especially in bn and b 33 ) in the outer-layer (y/8 > 0 . 75 ) in response 
to the relaxation in the shear rate. Figure 8 shows a comparison of the profiles of Sk/e obtained 
from the DNS results of Kim (unpublished), and predictions with the k-e model and the SSG 
Reynolds stress closure model. These all show that Sk/e is nearly constant in the inner layer and 
starts to decay rapidly beyond y/8 = 0.7. All the models show only a mild decay in the anisotropy 
level of the 1 1 and 33 components in this region, but DNS and experimental data presented in 
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Fig. 9 all show strong relaxation towards isotropy in correspondence with the decay of the shear 
rate. The main effects of increasing the Reynolds number are seen to be the reduction of the 
anisotropy close to the wall and a faster rate of return to isotropy near the center of the channel. 
Figure 7 shows that all the pressure-strain models (in conjunction with the MH diffusion model) 
produce reasonably good prediction of the turbulent kinetic energy. 

4. Concluding Remarks 

The k-e model and the second-moment closure models, apart from the SL model, produce 
similar predictions of the mean flow velocity, which agree well with the logarithmic law of the 
wall over most of the channel cross-section. In agreement with experimental observations in 
channel flows, there is only a small wake component. The SL model underpredicts the shear 
stress, and this leads to poor prediction of the mean flow velocity. All the other models predict 
the shear stress distribution correctly. 

The different models for the turbulent diffusion have little effect on the normal components 
of the Reynolds stress tensor in the log-layer, but strongly influence the rate of relaxation towards 
isotropy in the outer layer near the center of the channel. The MH diffusion model gives the 
best agreement with experimental data. 

LRR and SSG pressure- strain models give the best prediction of the streamwise and 
transverse components of the Reynolds stresses. The LRR model requires wall-reflection terms 
to achieve this but the SSG model does not. The wall reflection terms remain pronounced in 
the outer layer. Such an outer-layer influence is perhaps physically inappropriate. Only the SSG 
model could predict the lateral component of the Reynolds stress anisotropy correctly. 

The models failed to predict correctly the rate of relaxation of the streamwise and lateral 
components of the Reynolds stresses towards isotropy when the shear rate decreases in the outer 
layer y/6 > 0.7. 
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The pressure strain models which performed well in the present study, such as the LRR and 
SSG models, are expected to give reasonable results in other equilibrium shear flows. But there is 
no indication as to their performance in other complex flows. That is a subject for future studies. 
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Table 1 Diffusion terms Dy for developed thin shear flow 
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Figure 1 Geometry and coordinate system for plane channel flow. 
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Figure 2 Comparison of velocity profiles computed with the SSG pressure- 
strain model and various diffusion models to the log— law. 
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Figure 3 Comparison of components of the Reynolds stress anisotropy tensor 
computed with the SSG pressure-strain model and three diffusion models 
to experimental data of Laufer (1951). 
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Figure 4 Comparison of profiles of the turbulent energy computed with the 
SSG pressure-strain model and various diffusion models to experimental 
data of Laufer (1951). 
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Figure 5 Comparison of velocity profiles computed with the MH diffusion 
model and various pressure— strain models to the log— law. 
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Figure 6 Comparison of components of the Reynolds stress anisotropy tensor 
computed with the MH diffusion model and various pressure— strain 
models to experimental data of Laufer (1951). 
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Figure 7 Comparison of profiles of the turbulent energy computed with the 
MH diffusion model and various pressure -strain models to experimental 

data of Laufer (1951). 
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Figure 9 Comparison of Reynolds stress anisotropy tensor in channel flow 
from DNS of KMM (1987) and Kim (unpublished) and experiments of Laufer (1951 


27 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0 1 88 


5. FUNDING NUMBERS 

C NAS1-18605 


WU 505-90-52-01 


sj§i 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE | 3. REPORT TYPE AND DATES COVERED 

Apri 


4. TITLE AND SUBTITLE 

SYSTEMATIC STUDY OF REYNOLDS STRESS CLOSURE MODELS 
IN THE COMPUTATIONS OF PLANE CHANNEL FLOWS 


6. AUTHOR(S) 

A, 0. Demur en 
S. Sarkar 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665-5225 


9. SPONSORING /MONITORING AGENCY NAME{S) AND ADDRESS(ES) 

National Aeronautics and Space Adminstration 
Langley Research Center 
Hampton, VA 23665-5225 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No. 92-19 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA CR- 189646 
ICASE Report No. 92-19 


11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Michael F. Card 
Final Report 


12a. DISTRIBUTION/ AVAILABILITY STATEMENT 
Unclassified - Unlimited 


Submitted to ASME Journal 
of Fluids Engineering 


| 12b. DISTRIBUTION CODE 


Subject Category 34 
13. ABSTRACT (Maximum 200 words) 

This paper investigates the roles of pressure-strain and turbulent diffusion models 
in the numerical calculation of turbulent plane channel flows with second-moment 
closure models. Three turbulent diffusion and five pressure-strain models are util- 
ized in the computations. The main characteristics of the mean flow and the turbu- 
lent fields are compared against experimental data. All the features of the mean 
flow are correctly predicted by all but one of the Reynolds stress closure models. 

The Reynolds stress anisotropies in the log layer are predicted to varying degrees 
of accuracy (good to fair) by the models. None of the models could predict correctly 
the extent of relaxation towards isotropy in the wake region near the center of the 
channel. Results from the direct numerical simulation are used to further clarify 
this behaviour of the models. 


I 14. SUBJECT TERMS 
‘ turbulence modeling 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 


15. NUMBER OF PAGES 

29 


16. PRICE CODE 

A03 


| 20. LIMITATION OF ABSTRACT . 


Standard Form 298 (Rev 2-89) 

Prescribed by ANSI Std 239- <8 


NASA Langley. 1992 




